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Abstract 

The Bose-Einstein condensates recently created in trapped atomic 
gases are mesoscopic systems, in two senses: (a) Their size fall between 
macroscopic and microscopic systems; (b) They have a quantum phase 
that can be manipulated in experiments. We review the theoretical 
and experimental facts about trapped atomic gases, and give examples 
that emphasize their mesoscopic characters. One is the dynamics of 
collapse of a condensate with attractive interactions. The other is the 
creation of a ID kink soliton that can be used as a mode-locked atom 
laser. 



1 Mesoscopic systems 

A mesoscopic system has two distinctive characteristics: 

• It contains many more particles than a microscopic system, but is small 
enough for us to manipulate in the laboratory. 

• It possesses a quantum phase that we can control. 

The term "mesoscopic" has been associated with "quantum dots," in 
which the passage of electrons can be controlled singly. Now we have new 

*Expanded version of an invited talk given at the third joint meeting of Chinese physi- 
cists worldwide, Chinese University of Hong Kong, August, 2000. 
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mesoscopic systems, with the achievement of Bose-Einstein condensation in 
a system of trapped atoms. These Bose condensates typically contain the 
order of 10 6 atoms, in a magnetic or optical trap, at a temperature of the 
order of 10 _7 K. After a brief summary of known facts, we shall give examples 
of both aspects mentioned above. For more details and literature, we refer 
to recent review papers on experiments and theory [Q. 

2 The Bose condensate 

As we know, there are two types of elementary particles: fermions of half- 
integer spin, which cannot occupy the same state, and boson of integer spin, 
which love to be in the same state. Bound states of elementary particles, 
such as atoms and molecules, also fall into these categories. The bosonic 
atoms we are dealing with include 87 Rb, 23 Na, 7 Li, and atomic hydrogen. 

In the ground states all bosonic atoms occupy the same state; forming 
the Bose-Einstein condensate. What makes this an identifiable entity is its 
stability, arising from the quantum-mechanical fact that atoms are indistin- 
guishable. In Fig. [I] we illustrate how this reduces the number of excited 
states. The situations illustrated would represent three different states, if 
the particles were distinguishable. But for indistinguishable particles they 
are one and same the same state. Because of this, there are fewer ways to 
excite the condensate, and hence it has a high degree of stability. 

Ground state: 

ooooooooooo 

ExL-ilyd slates: 

o°ooooooooo 

000*0000000 

ooooooooooo 

Classical: 3 different states 
Quantum: Same state 

Figure 1: The indistinguishability of particles reduces the number of excited 
states, and makes the condensate stable. 

In fact, when the temperature is raised from absolute zero, there remains 
a finite fraction of particles in the condensate. The condensate fraction de- 
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creases with rising temperature, and disappears only above a finite temper- 
ature T c , the critical temperature of the Bose-Einstein condensation. 

In contrast, for "Boltzmann statistics", the ground state wave function is 
the same, but there will be a high density of excited states. In this case, no 
condensate exists at any finite temperature. That is, the transition temper- 
ature is at absolute zero. 

The process of Bose-Einstein condensation is illustrated in Fig. ^| 
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Figure 2: The approach to Bose-Einstein condensation. 



3 Experimental techniques 

The process to achieve BEC in the laboratory is illustrated in Fig. An 
atomic beam was slowed down by laser cooling, and a magnetic trap was 
turned on. The spin up atoms are confined, while the spin down one are 
expelled. Final cooling is achieved through evaporation and rethermalization. 
To evaporate the fast atoms, rf light of appropriate frequency is turned on 
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Figure 3: Experimental steps in the trapping and cooling of atoms. 



to induce spin flip. In the final condensate, we have typically 

iV~ 10 6 

V ~ KT 6 cm 3 

- ~ 10 12 cm" 3 (I) 

The density is some seven orders of magnitude smaller than that of air in 
the atmosphere. 

A brief history of Bose-Einstein condensation in trapped gases is given 
Fig-S 



4 Helium clusters 

Long before trapped atoms, of course, Bose-Einstein condensation had been 
seen, though indirectly, in the superfluidity of liquid 4 He. Andronikashvili 
performed the classic experiment on the superfluid fraction, by measuring the 
moment of inertia of a stack of rotating plates immersed in liquid helium, 
as a function of temperature. The idea was that the normal fluid would be 
dragged into rotation with the plates, but the superfluid component would 
not participate in the motion. 
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History 
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Figure 4: History of atom trapping and cooling. 

The ultimate miniaturization of the Andronikashvili experiment has now 
been done in 4 He clusters of a few hundred atoms, with the rotating plates 
replaced by a rotating atom of 3 He. Superfluidity was reported in a cluster 
with 60 atoms. This is a very interesting mesoscopic system, but beyond the 
scope of our discussion here. 

5 Quantum coherence 

For an ideal Bose gas, the condensate is made up of particles in the same 
single-particle quantum state. When there are interparticle interactions, 
single-particle states are no longer meaningful; but we can define the conden- 
sate wave function as the quantum amplitude for removing a particle from 
the ground state: 

V>(r, t) = yfce* (2) 

where n s is the condensate density, and <fi is the quantum phase of the con- 
densate. 

The quantum coherence of the condensate is a remarkable thing. In 
experimental traps, the condensate can have a spatial extension as large as 1 
mm, which is almost visible to the naked eye, and yet this glob of matter is 
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characterized by a single wave function, with a single phase. This has been 
demonstrated experimentally by switching off the magnetic trap, cutting 
the condensate in half with a laser "knife," and allowing the two halves to 
fall under gravity. The two condensates expand as they drop, and eventual 
overlap in space, creating a interference pattern, as illustrated in Fig. [5]. An 
actual photograph of the interference fringes in shown in Fig. |6|. 




Figure 5: A condensate in free fall is cut by a laser knife, and the two halves 
continue to fall, expand, and overlap. 




h'ringe .spacing = k ! 2 = 1 5jum 



Figure 6: Interference fringes in two overlapping condensates of sodium 
atoms. {Courtesy W. Ketterle.) 

We can describe two non-overlapping free condensates by wave functions 
that are nearly plane waves: 

fa = Cie i(^-r-"it) 

fa = c 2 e i{ ^ 2 - r - W2t) (3) 
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When they overlap, they become one condensate with wave function 

= V»i + ^2 (4) 
The density of the system is given by 

|^| 2 = |d| 2 + |Ci| 2 - 2Re [C*C 2e i(ki-k 2 )-r-i(a; 1 -^)t)j ^ 

where the last term exhibits the interference fringes. From Fig. j|, the fringe 
spacing gives a half wavelength A/2 =15 /im, which corresponds to a relative 
wave number 

h - h = y (6) 
The relative velocity of the two condensates was therefore given by 

h 

v\ — V2 = — (ki — k 2 ) = 0.05 cm/s (7) 
m 

Atoms in the condensate in a magnetic trap can be ejected by inducing 
spin flips, through application of rf light. The expelled atoms fall under 
gravity as pulses of coherent matter, as shown in Fig. [7[ This is the first 
example of an atom laser ||. 




Figure 7: The first atom laser. {Courtesy W. Ketterle.) 
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6 Nonlinear Schrodinger equation 



The condensate wave function satisfies a nonlinear Schrodinger equation 
(NLSE), also called the Gross-Pitaveskii equation, in a mean-field approxi- 
mation that neglects the coupling to the thermal cloud: 

i *^=(-^ + v<r) + 9\1>\ 2 )l> (8) 



with 

A^U 2 a , nS 

9 = 9 

m 

where a is the S-wave scattering length. The normalization condition 

d 3 r\ij\ 2 = N (10) 



determines the number of particles in the condensate. This is a constant of 
the motion, but its value cannot be changed by convention, contrary to that 
for an ordinary wave function, because the equation is nonlinear. 



7 The scattering length 

The scattering length a is illustrated in Fig. |8|. The reduced relative wave 
function of two atoms, as a function of the separation r, approaches a linear 
asymptote beyond the range of the potential. The scattering length is the 
the intercept of the asymptote on the r axis. 




U — Acaciwa pocrtra 



Figure 8: Scattering length a for potentials without bound states. 
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When the potential has no bound state, as illustrated in Fig. [| a repulsive 
potential give a > 0, and an attractive potential give a < 0. When there 
are bound states in a attractive potential, however, the scattering length 
can have either sign. As the potential gets more attractive, it changes sign 
whenever a bound state occurs. 




J 

Figure 9: Scattering length a for potentials without bound states. 

One of the exciting recent developments is that the scattering length can 
be tuned, from — oo to oo, by changing an external magnetic field, through 
the Feshbach resonance in atomic scattering, which is similar to a bound 
state 0. This realizes the usual thought experiment of changing the coupling 
constant in the Hamiltonian. 

It should be emphasized that a is a parameter that describes scattering in 
3D in the limit of zero energy, where the shape of the potential is irrelevant. 
At higher energies, the shape of the potential does make a difference, and 
one must include other parameters, such as the effective range 0. 



8 Quantum tunneling 

The trapped atoms that form a Bose condensate at such low temperatures 
would have solidified in free space. In the trap, the zero-point motion kept 
them apart. In the case of 6 Li, which has negative scattering length, the 
condensate is metastable, and decays via quantum tunneling, if the total 
number of atoms N is sufficiently small. If N exceeded a critical number, 
however, the system collapses @, || . Let us discuss each of these scenarios. 
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For negative scattering length, we have g < 0. Let us denote 

G = -g- 



Anh 2 a 



m 



(11) 



We can see qualitatively that in free space the condensate collapses in spatial 
dimension D > 2. The energy density of system is 

£ = !hwI 2 -gH 4 (12) 

2m 

Suppose the system is confined to a region of dimension R. Then ip\ 2 ~ 
N/R D , so 



E 



h 2 N 



R 2D 



(13) 



2m R D + 2 

Thus E -oo when R -> 0, if D > 2. 

If there is a confining potential, then collapse happens only if iV is suf- 
ficiently large, because there is a potential barrier against collapse. To see 
this, consider the Hamiltonian corresponding to (B: 

" *2 p -1 

^*V 2 V> + V{r)ip*ip {ip*ip) 2 (14) 



H\ 



2m 



and choose a harmonic oscillator potential 



V{r) 



x 2 2 

-mu r 



(15) 



For qualitative arguments, let us consider a confined wave function of a 
Gaussian shape 



ipo(r) = C exp (-ar 2 /<io) 



where 



do 



muj 



(16) 



(17) 



The parameter a describes the width of the Gaussian. A Gaussian narrower 
than the harmonic oscillator wave function corresponds to a > 1. We can 
then replace — |^V 2 ?/> — > x( r )V ; o, where 



X{r) 



2md\ 



3a — 



a 



(18) 
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Thus the Hamiltonian can be represented in the form 



HI 



d 3 rtt r (ip) 



(19) 



where the "field potential" £l r (i[)) is given by 
n r ^) = [V(r) + x(r)] |V 



2 ^l./.|4 



(20) 



This potential is plotted as a function of for different rT, in the right 
panel of Fig. [TO. There is an energy barrier with height 



W r = ^[V(r) + X (r)] 2 



(21) 



which is nonzero at r = 0, because x(0) 7^ 0. At large \ip\ the field potential 
tends to —00. In reality, of course, the NLSE ceases to be valid somewhere 
along the drop, for other physical effects, such as solidification, come in. 



Condensate 
wave function 



Field potential 




r > 



Figure 10: The field potential as seen by the condensate wave function. 
Local quantum tunneling A — > B occurs across an energy barrier. If the 
wave function at r goes over the energy barrier, then there is local collapse 
at that point. 



A typical initial wave function ipo(r) is shown on the left panel of Fig. TH. 
At a given r, we can measure the wave function on the left panel, and transfer 
it to the horizontal axis on the right panel. If it lands to the left of the barrier 
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maximum for the particular r, then the system at that point is classically 
stable, but can decay via quantum tunneling as indicated by the classically 
forbidden path A — > B. If it lands to the right, the system at that r will 
rapidly collapse to a state of high density. 

Since ipo(r) has a maximum at r = with value proportional to yN, 
the system will decays local tunneling, if it is classically stable at r = 0. 
Otherwise, the system in a region about r = will rapidly collapse. The 
condition for stability against collapse is therefore f2oC0o(O)) < Wo. Within 
the Gaussian assumption, this corresponds to iV < iV c , with 

Nc = 3 JWdo = 0M5 p_ 
8 V ct \a\ y/a \a\ 

In the experiments with 7 Li at Rice University |§, the parameters are 

a = —1.45 nm 
do = 3 m 

2tt , , 

= 6000s 23 

OJ 



With the choice a = 1.8 (See later), we obtain N c ~ 1200 from 

For N < N c , the condensate decays slowly via quantum tunneling, as 
indicate by the path A — > B in Fig. [II]. Since the tunneling path depends 
on r, the tunneling is a local phenomenon, proceeding at different rates at 
different distances from the center of the trap. 

It should be noted that quantum tunneling is a phenomenon not included 
in the NLSE, which is a classical approximation. The tunneling probability 
can be calculated via a Feynman path integral, and in a WKB-type approx- 
imation can be obtained by solving the NLSE in imaginary time 0. 

Qualitatively, we can describe the phenomenon as follows. The conden- 
sate locally makes a quantum transition from a low-density state A to a 
high-density state B. The time at which this happens is random, and follows 
a statistical distribution. Since the total number of particles is conserved, 
what we see is that high-density droplets would "ooze out" of the conden- 



sate, as illustrated in Fig. 11. Once form, the droplets will tend to contract 
to states of even higher density, and will leave the trap through effects not 
included in our Hamiltonian, such as spin-flip scattering. 

A more rigorous treatment of the field potential, independent of the Gaus- 
sian assumption, may be modeled after that of the "effective potential" in 
quantum field theory [|K|. 
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Figure 11: High-density droplets "ooze out" of the condensate via local quan- 
tum tunneling. 

It should be noted that, unless iV « N c , the tunneling probability is very 
small. In reality, it would be masked by the decay of the condensate through 
spin-flip scattering, an effect that we have not taken into account. 



9 Collapse of the condensate 

When N > N c , the condensate with negative scattering length collapses. 



As we can see from Fig. [1CJ, the collapse does not happen all at once over 
the entire condensate, but begins at the center of the condensate, where the 
density is highest. A black hole opens up at the center, in which, theoretically 
according to the NLSE, the density fluctuations becomes infinite, and the 
pressure becomes negative |7|]. 

Numerical solutions of the NLSE have been obtained f|, with the perime- 
ters (^). After the condensate was created in some arbitrary state, the wave 
function quickly adjusted itself to a Gaussian-like quasi-stationary state, with 
a width corresponding to a = 1.8. This situation persists infinitely if N is 
sufficiently small. 

If N is larger than a critical value, then the quasi-stationary states lasts 
for a time, but then suddenly a finite volume about the center implodes, 
creating a black hole. Thereafter, waves of implosion arrives from all over 
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Figure 12: The condensate is stationary for a time, and suddenly collapses 
locally, near the center of the trap, opening up a black hole. The density in 
the black hole shoots up as time goes on, fed by waves of implosion, which 
make the density look turbulent. The density height is truncated because of 
finite grid size in the numerical simulation. 

the condensate, whose density now appears to be turbulent. These features 
are summarized in Fig. [12], where r\ip\ is plotted as a function of r in unit of 
do, and t in unit of 2/u, for a fixed N. By decreasing N until the time of 
collapse goes to infinity, it was found that N c rs 1200. 

In the experimental situation, the condensate in the trap can exchange 
atoms with an uncondensed thermal cloud, which contains far more atoms 
than the condensate. Equilibrium is reached when the chemical potentials 
equalize. The atoms can be depleted through two and three-body spin-flip 
collisions, and the latter becomes especially important where the density 
becomes large. These gain and loss mechanisms can be simulated by adding 
extra pure-imaginary terms to the NLSE, with a linear gain term, and a loss 
term proportional to |?/'| 4 V ; - With these additions, the system settles down 
to steady state in which N oscillates in growth-collapse cycles, as shown in 
Fig. [L3| Again, we obtain N c « 1200. 
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Time (s) 

Figure 13: The condensate goes through steady-state growth-collapse cycles, 
in the presence of gain-loss mechanisms. 

10 Solitons in a condensate ring 



A spatial variation of the phase <fi of the condensate wave function gives rise 
to a superffuid velocity 



and a super current density 



v s = — V0 

m 



3s = n s v s 



(24) 



(25) 



Thus, flow patterns in the condensate can be created by "phase engineering." 
This has been done to create vortices jXTJ, |T2] and solitons |TB| . 



One-dimensional solitons have been used in optical fibers for signal trans- 
mission over thousands of kilometers, with minimal loss. The electromagnetic 
field in an optical fiber obeys a NLSE of the same form as (ffl) Thus, 



it is natural to wonder whether we can create solitons in a ID condensate. 
As illustrated in Fig. Ill, the nonlinear term g\i/j\ 2 represents an effective po- 



tential, which is attractive for g < 0, and repulsive for g > 0. g < OAs we 
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saw earlier, the case does not lead to collapse in ID. In this localized 
wave packet created a attractive well of the same shape, hence it would trap 
itself into a bright soliton. In the case g > 0, on the other hand, a "kink" 
configuration would self-trap, giving rise to a dark soliton. 




Figure 14: Formation of ID solitons via "self-trapping". 



Consider now a ring-shaped condensate. If the radius of the ring is much 
larger than the radius of the cross section, then we can regard the system 
as one-dimensional, in the sense that transverse excitations have a much 
high energy than the logitudinal ones, and may be ignored at low tempera- 
tures. However, the cross sectional radius may still be large compared to the 
scattering length, so that the system is still 3D with respect to interatomic 
scattering. 

We concentrate on the repulsive case with g > 0. The solution to the 
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NLSE must satisfy the boundary condition 



ij(9 + 2n,t) =ifi(6,t) (26) 

where 9 is the angle around the ring. Let us write 

V(M =mt)e^ (27) 

For a kink soliton, / must have a zero, and so changes sign when 9 increases 
by 27r. To satisfy the boundary condition, the phase must change by an odd 
multiple of tc in one complete revolution. The simplest case corresponds to 
a change of 7r in the phase. This means that the lowest kink soliton has 
angular momentum per particle h/2, for the same mathematical reason that 
an electron has spin 1/2. 



The solution can be obtained analytically [16]. The phase changes most 
rapidly at the location of the kink, and its slope there gives the kink a nonzero 
propagation velocity. Thus, the soliton is manifested as a dark spot moving 
around the ring with a characteristic velocity. The width of the soliton turns 
out to be proportional to iV" 1 / 2 , where N is the number of atoms in the 
condensate. 



11 A mode- locked atom laser 

The first atom laser illustrated in Fig. [7| consists of pulses of coherent atoms. 
However, the different pulses do not have definite phase with respect to each 
other. There is an advantage to have a mode-locked laser — one in which 
the output pulses are phase coherent. In more recent experiments, a steady 



stream of output pulses has been achieved |jl5| , but it is not clear whether 
the pulses are coherent with respect to each other. 

We can use of our kink soliton to make a mode-locked atom laser Jl6| 



The physics is illustrated in Fig. [TR The ring of condensate is surround by 



a thermal cloud, which tends to maintain an equilibrium number of particles 
N in the ring. Now suppose we use some stroboscopic device to take particles 
out the soliton, by enlarging its width, which is proportional to iV -1 / 2 . Then 
particles will flow from the thermal cloud to the condensate to restore N, 
and thus restoring the width. 

An important question is how to create the kink soliton, and whether it 
is robust against perturbations. In a simulation, the soliton was created by 
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Figure 15: Mode- locked atom laser created by stroboscopically shaving off 
atoms at the kink of the dark soliton. The width grows back due to gain 
from the thermal cloud. 



applying a moving momentary repulsive potential, such as that created by a 
laser beam, to create a hole in the condensate that moves at the character- 
istic soliton speed. This hole develops into the kink soliton, plus all kinds of 
excitations such as sound waves. To clean up the configuration, a periodic 
loss mechanism was introduced, in the form of am imaginary stroboscopic 
potential that also acts as outcoupler. At the same time, a gain mechanism 
was put in to simulate the thermal cloud. All unwanted excitations rapidly 
disappeared, leaving the kink soliton in steady-state equilibrium. The sta- 
bility of the soliton come from topology — the fact that it has half-integer 
angular momentum per particle. 

The creation and stabilization of the soliton is depicted in Fig. [16], where 
\ip | is plotted as a function of 6 and t, with a contour plot projected onto the 
9-t plane. The fact that we have in fact a kink soliton is indicted in Fig. 17] , 
where the angular momentum per particle is plotted as a function of time. 
Finally, Fig. |TJ^ shows the oscillation of the number of particles in the ring. 
The particles lost from the ring form the pulses in the mode-locked laser. 
They are coherent because the stroboscopic action is coherent. 
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Figure 16: The kink soliton was created by stirring the condensate at its 
characteristic velocity, Other excitations also created in the stirring were 
cleaned up by a stroboscopic loss mechanism, which also serves as output 
coupler for the atom laser. The kink soliton survives because of topological 
stability. 
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